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ABSTRACT 

We study star formation history and stellar metallicity distribution in galaxies in a A cold dark matter 
universe using a hydrodynamic cosmological simulation. Our model predicts star formation rate declining 
in time exponentially from an early epoch to the present with the time-scale of 6 Gyr, which is consistent 
with the empirical Madau plot with modest dust obscuration. Star formation in L* galaxies continues 
intermittently to the present also with an exponentially declining rate of a similar time-scale, whereas 
in small galaxies star formation ceases at an early epoch. The mean age of the extant stars decreases 
only slowly with increasing redshift, and exceeds 1 Gyr at z — 3. Normal galaxies contain stars with a 
wide range of metallicity and age: stars formed at z < 1 have metallicity of 0.1 — l.0Z@, while old stars 
take a wide range of values from 10~ 6 Zq to 3.0Zq. The mean metallicity of normal galaxies is in the 
range 0.1 — I.OZq. Dwarf galaxies that contain only old stars have a wide range of mean metallicity 
(10 -4 — I.QZq), but on average they are metal deficient compared with normal galaxies. 

Subject headings: stars: formation — galaxies: formation — cosmology: theory — methods: numerical 



1. INTRODUCTION 



Over the last 15 years, cold dark matter (CDM) mod- 
els have served as basic frameworks to study the forma- 
tion of cosmic structure, and have been successful in delin- 



eating a scenario of gala xy formation ( Davis et al., 1985 
Blumenthal et al., 1984). We now understand the forma- 



tion of large-scale structure reasonably well in terms of 



(Efstathiou, Sutherland, & Maddox, 1990; 


Ostriker anc 


Stcinhardt, 19951; Turner and White, 1997 


iBalbi et al. 


2000; Lange et al., 200q; Hu, et al., 2000). 



When and how galaxies formed, however, are a far 
more complicated problem due to non-gravitational phys- 
ical processes operating on small scales. Many authors 
use semianalytic models to study galaxy formation, where 
dark matter halo formation under the hierarchical struc- 
ture formation is supp lemented hcuristically with physical 



wc obtain meaningful results on some aspects of galaxy 
properties. It is certainly impossible to resolve internal 
structure of galaxies with the present simulations. Nev- 
ertheless, we would expect that global properties such as 
global star formation rate (SFR) are described reasonably 
well with a few adjustable parameters, since they primar- 
ily depend on thermal balance of the bulk of clouds, such 
as how gas cools, or how surrounding gas is reheated by 
feedback by star formation. The fact that we obtain a rea- 
sonable amount of baryonic mass frozen into stars and a 
global metal abundance with reasonable input parameters 
justifies, at least in part, our expectation. 

In this paper, we primarily discuss the star formation 
history of galaxies and the evolution of stellar metallicity, 
which are the direct output of the simulation. The aspects 
that require additional use of a stellar population synthe- 
sis model, such as luminosity function and colo rs of galax- 



1998 



processe s for baryons ( White fc Frcnk, £ Daugh ct -grl 



Kauffmann ot al., 1999a 



ies, will be discussed in a separate p ublication Nagamine 



1999| ; |Cok~ et al., 2000| ). They have succeeded to give a 
picture of galaxy formation roughly consistent with obser- 
vations, but have the disadvantage that one has to assume 
a set of simplified model equations, often associated with 
additional free parameters, for each physical process in- 
cluded in the model. The advantage of semianalytic mod- 
els is, on the other hand, that they are computationally 
light, and therefore, one can search for a viable model with 
a small amount of cost. 

An alternative but computationally much more expen- 
sive approach is direct cosmological hydrodynamic simu- 
lations. So far, intergalactic medium and overall galaxy 
formation processes h ave been well studied wi t h hydrody- 
namic approach (e.g., |Ccn and Ostriker, 1992a , 1993 , 2000 



h ! HI j A ■ T"l Fukugita, Cen. and Ostriker (2001, hereafter Paper II 
homorvillo and Primack, " ° . '—. — , i — r -?' r 

I -. — : : ' Smrn rnp hvnrnnvnflTnir flr-mrnnr'h i« pnmmirfirinnn v ni 



Katz, Weinberg, fe Hcrnquist, 1996] ) (sec |Pcarce et al 



1999, however). Much effort invested in improving the ac- 



curacy of simulations has brought the hydrodynamic mesh 
approaching to (1000) 3 , with which we may hope that 



Since the hydrodynamic approach is computationally de- 
manding, we do not attempt to make a fine tuning of input 
parameters, but our aim here instead is to make qualitative 
predictions which can be used to either verify or falsify the 
ACDM model by future observations with only a limited 
number of assumptions, treating the dynamics of baryons 
more accurately than in semianalytic models. 

Our results would give insight on how galaxies are as- 
sembled from an early epoch to the present. Observation- 
ally interesting questions would include: (i) How old are 
stars in normal and dwarf galaxies, and what is their age 
distribution? A relevant question is whether star forma- 
tion takes place continuously or intermittently; (ii) What 
is the stellar metallicity distribution in normal and dwarf 
galaxies? A relevant old problem is the paucity of metal 
poor stars in the solar neighborho od (G-dwarf problem: 
van den Bcrgh, 1962 ; Schmidt, 1963 ); (iii) Is there a unique 
age-metallicity relation for stars in normal galaxies, and 
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(iv) Is there a relation between the mass of galaxies and av- 
erage metallicity? How these aspects evolve with redshift 
is also an observationally relevant problem, although the 
observational studies could answer only for global prop- 
erties for high-redshift objects. Our simulation yield un- 
ambiguous predictions to these problems, at least qualita- 
tively. 

In § ||, we describe our simulation. In § ^ we first present 
the global star formation history of the entire simulation 
box. We then discuss in § || the star formation history 
of individual galaxies. The mean age of the stars and 
the formation epoch of galaxies are discussed in § ||. We 
discu ss the stellar metallicity distribution and the mean 



metal icity of galaxies in § 6 , and conclude in § 7[ 



2. SIMULATION AND PARAMETERS 

We use a recent Eulerian hydrodynamic cosmological 
simulation with a comoving box size of 25ft. -1 Mpc with 
768 3 grid cells and 384 3 dark matter particles of mass 
2.03 x lO 7 /^ 1 M Q . The comoving cell size is 32.6/i _1 
kpc, and the mean haryonic mass per cell is 3.35 x 
10 5 /i~ | 1 M<t i (h is the Hubble constant in units of Hn = 



100 km s 1 Mpc 1 ). The cosmological parameters are cho- 
sen to be O m = 0.3, n A = 0.7, VL b h 2 = 0.016, h = 0.67, 
as = 0.9, and the spectral index of the primordial mass 
power spectrum n — 1.0. The present age of the universe 
is 14.1 Gyr. The structure of t he code is similar to that of 



Cen and Ostriker (1992a , 1995 ) , but significantly improved 



over the years 

The Eulerian scheme we use here has higher mass res- 
olution than other hydrodynamic approaches, su ch as 



the smoothed particle hydrodynamics (SPH; e.g., Katz 
Weinberg^ fc Hcrnquist, 1996] ; Dave, Dubinski, & Hcrn- 
quist, 1997 | ) and the adaptive mesh refinement (AMR; e.g 



Bryan fc Norman, 1995 ; Kravtsov, Klypin, fc Khokhlov, 
1997[ ), while SPH and AMR usually have higher spatial 



resolution than the Eulerian method. We remark that our 
code is specially designed to capture the shock region well, 
using the total variation diminishing method ( Ryu et al.J 
1993). The higher mass resolution means more resolving 
p ower, at higher wavenumbers for a given power spectrum. 



and hknce a, better treatment of early structure formation. 

The baryons are turned into stars in a cell of overden- 
sity Sp/p > 5.5 once the following three conditions are 
satisfied: (i) contracting flow V • v < 0, (ii) fast cooling 
icooi < £dyn, and (iii) Jeans instability TO gas > mj. Each 
stellar particle has a number of attributes at birth, includ- 
ing position, velocity, formation time, mass, and metallic- 
ity. Upon its formation, the mass of a stellar particle is 
determined by m t — c x m v ,,yA.tJt^ where At is the cv 



galaxy. The stellar particle is placed at the center of the 
cell after its formation with a velocity equal to the mean 
velocity of the gas, and followed by the particle-mesh code 
thereafter as collisionless particles gravitationally coupled 
with dark matter and gas. 

Feedback processes such as ultra-violet (UV) ionizing 
field, supernova (SN) energy, and metal ejection are also 
included self-consistently. The SN and the UV feedback 
from young stars are treated as follows: AE?sn = m»c 2 esN 
and AEuv = m»c 2 euvffiy, where g v is the normalized 
spect rum of a yo ung, Orion-like stellar association taken 
from |5calo (1986 ) , and the efficiency p arameters are taken 
as (e sN,euv) — (10~ 5 ,3 x 10~ 6 ) (cf. Gen and Ostriker, 
1993). The A_Esn is added locally in a cell, and the 



Ai?uv is averaged over the box. Extra UV is coadded 
for quasars as A£agn = m*c 2 eAGN.fi^ where f u is an 



AGN spectrum taken from Edelson and Malkan (1986), 
and eagn = 5x 10~ 6 is adopted. The proportionality be- 
tween Ai/AGN and m* is an assumption, and our ignorance 
is included in eagn which is adjusted to fit the hard X-ray 
backgroun d observed by ASCA (see Phillips, Ostriker, & 
Cen, 2000). Metals are created according to mz = Ym*, 
where mz is the mass of metals, and Y is the yield. The 
fraction 25% (f e j = 0.25) of the initially collapsed baryons 
is ejected back into the intergalactic medium in the local 
grid cell. This ejected g as is polluted by metals, and 2% 
(Y = 0.02; Arnett, 1996) of initially collapsed baryons are 
returned to intergalactic medium as metals, which are fol- 
lowed as a separate variable by the same hydro-code which 
follows the gas density. The value of the three parameters 
(csn, £uv, Y) are proportional to the amount of high-mass 
stars for a given amount of collapsed baryons. The results 
presented in this paper is insensitive to the assumed value 
of f e j provided that it is non zero. Thus, all our ignorance 
concerning star formation is parameterized into essentially 
one free parameter, the efficiency of cooling, i.e., collaps- 
ing matter that is transformed into high-mass stars, and 
this value is empirically calibrated so that one obtains a 
final cluster gas-metallicity of Z = hZ Q , where Z Q — 0.02 
is the solar metallicity. 

We use the HOP grouping algorithm ( Eiscnstcin and| 



Hut, 1998) to identify galaxies with threshold parameters 
(Pouter, Saddle, (5 P cak) = (80,200,240). We set the mini- 
mum number of grouping particles to 5; changing this to 
2 does not modify the catalog. We identify 2097 galaxies 
in the entire box at z = 0, the minimum galaxy stellar 
mass being lO 67 /)." 1 M©. The total mass of stars corre- 
sponds to = 0.0052, i.e., 15% of baryons are collapsed 
into stars (including the ejected gas). This is consistent 



rent tlime-step in the simulation and m is the hary- 



with the upper limit of an e mpirical estimate of Fukugita, 
Hogan, and Peebles (199^ ). The average metal density 



onic gas mass in the cell. We take t» = max(£dyn, 10 7 yrs) 
where 10 7 yrs is the shortest time-scale of star formation 
for O stars. We assume the star formation efficiency to be 
c* = 0.25. The fraction of baryonic mass that collapses 
into dense state is fairly well defined by our code, nearly 
independent of the detailed values of c» and the minimum 
value of t*, as confirmed by a recent numerical experiment 
of Norman (2000, private communication). 

The mass of the stellar particles ranges from 10 35 — 
lO 103 /^ 1 M©; with the mean (m*) = lO 69 /^ 1 M© at 
z = 0. An aggregation of these particles is regarded as a 



is 10 71 M©Mpc , which is also close to the empirical 
estimate of 10 74 M© Mpc -3 (Fukugita & Peebles, unpub- 
lished), justifying our choice of parameters. 

We examine the merger history of galaxies, and find that 
93% of the galaxies at z = 0.5 preserve more than 80% of 
their stellar mass to z = 0: i.e., most of the constituent 
stellar particles are not lost upon merging or tidal shred- 
ding. This tells us that the galaxies identified in our simu- 
lation are dynamically stable enough to obtain meaningful 
results for the evolution of individual galaxies. 



Table 1 



Parameters of Exponential Fit for SFR 



Redshift 



Mass range [h 1 Mq] 


A [h 2 Mqyt- 1 Mpc~ 3 ] 


r [Gyr] 


All 


0.24 


5.9 


M atellar > 2 x 10 9 


0.23 


6.1 


2 X 10 8 < A^ellar < 2 X 10 9 


0.021 


0.83 


Mrtellar < 2 X 10 8 


0.011 


0.55 



NOTES. — Parameters ghjen above are the least-square-fit of the 
SFR in Figure 111 and to the following exponential function: 
SFR [h 2 M yr-*Mpc- 3 r= A exp(-t/r), where A is the SFR at 
t = 1 Gyr, and t gives the e-folding time of star formation. The rise 
in the first 1 Gyr is neglected upon fitting. 



3. GLOBAL STAR FORMATION RATE 

Wc first present in Figure la the global SFR as a func- 
tion of time, going from right to left. The prediction of 
global SF R in the CDM mode l has been discussed in the 
litera ture ( Baugh ct al., 1998 1_ Nagaminc, C en, and Os 



triker j 2000| ; Komcrvillc, Primack, fc Fabcr, 2001| ) \'\. This 
quantity is observationally known up to uncertainties as- 
sociated with dust obscuration, and the agreement with 
observations can be used as a verification of the calcu- 
lation. This quantity does not receive the effect of the 
grouping procedure, but gives a measure of the efficiency 
of gas cooling, which depends on the thermal balance. 

We also indicate a boxcar smoothed histogram with 10 
point running average with the dotted curve. The dashed 
line is a least-square fit to the exponential decay function, 
corresponding to a time-scale of t = 5.9 Gyr. 

The identical information is presented in Figure lb, but 
using the logarithmic scale for SFR and redshift as a time- 
scale in accordance with the Madau plot. In this diagram 
the meaning of solid and dotted curves are reversed: the 
dotted curve follow the raw data, and the solid line is 
the smoothed data. The observation is corrected for dust 



extinction according to the prescription of Stcidcl ct al 
(1999§: we assume highly uncertain extinction correction 
factors to be 1.3 [z < 2) and 2.5 [z > 2) to obtain a bet- 
ter fit, while Stcidcl ct al. (191)91 ) used higher values 2.7 
(z < 2) and 4.7 (z > 2). With a rather large error in 
the empirical estimate of SFR, a global agreement is seen 
between our calculation and the observation. We do not 
observe a peak of SFR at low redshift, as often found by 



semi analytic modellers of galaxy formatio n (Baugh et al. 
1998 fc jSomcrvillc, Primack, fe l abor. 200 l|i P. The SFR is 
on average a smooth function of redshift, and nearly levels 
off at z > 4. Note that this behavior is consistent with an 
exponential decay of the SFR in time. This means that 
the bulk of stars have formed at redshift higher than 2: 
25% of stars formed at z > 3.6, and another 25% were 



1 The simulation used by Nagamine et al. (2000) assumed a too 
high yield of metals (Y = O.Orj), which resulted in overproduction 
of stars compared to the empirical estimate. Also, the simulation 
mesh was coarser and one may suspect the effect of low resolution 

ult presented in this paper 



at above z ss 4. We frin-^rW that +ho ro 



Nagaminc et al. (2000) 



supersedes that of 

2 In semianalytic models, one may suppress star formation at high- 
redshift by taking a larger SN energy feedback parameters. In gen- 
eral, they require a strong SN feedback t.n fit th e faint-end of the 
galaxy luminosity funcition (see Cole et al., 2000). See Paper II for 



the luminosity function in our simulation. 
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Fig. 1. — (a) The SFR of the entire simulation box as a function of 
time, going from right to left. The dashed line is the least-square fit 
to the decaying exponential function given in Table h with a time- 
scale of t — 5.9 Gyr. (b) The identical information, the SFR, as a 
function of redshift (the Madau plot) The source of data points 
are given in Nagamine et al. (200C). Here we adopted factors of 
1.3 [z < 2) and 2.5 (z > 2) tor dust extinction correction of the 
observed data points. See text for discussion. 



added between z — 3.6 and 1.8 (By z = 1, 68% of stars 
were formed). We consider that our result is reliable up to 
z = 5 — 6 from a mesh effect consideration, whereas that at 
a yet higher redshift may receive an effect of poor resolu- 
tions. Observationally it is an unsettled problem w hether 
SFR levels off at high-red shift ([Steidcl ct al., 1999|) or still 
increases beyond z > 3 (Lanzetta et al., 1999). It is an 
interesting observational problem to find out how SFR be- 



4 



haves at high-redshifts, but this can be done only after 
proper knowledge of dust extinction or with the observa- 
tions that do not receive effects of extinction. For now, it 
is sufficient for us to know that the calculation of global 
SFR is grossly consistent with the observation. 

4. STAR FORMATION HISTORY IN GALAXIES 

In Figure || we present the star formation histories of 
galaxies in the simulation divided into three samples ac- 
cording to their stellar mass at z — 0. Note that the star- 
bursts at high-redshift may have taken place in smaller 
progenitors that later merged into a massive object that 
is classified in the largest mass bin. All galaxies which 
fall into each mass interval are co-added in the histogram. 
Most of the mass reside in large galaxies which are rep- 
resented in the top panel. This histogram can also be in- 
terpreted as the age distribution of stars in galaxies. The 
dashed lines show the least-square fits to the exponential 
decay function with parameters given in Table 0. The 
time-scale t = 6.1 Gyr for 'normal galaxies' (stellar mass 
Afgteiiar > 2 x 10 9 /i _1 M ; top panel of Figure 2) agrees 
qualitatively with the star formation history inferred for 
present-day Sbc galaxies, which are about the median type 
of local galax ies (Bandage, 1986; Gallagher, Hunter, & Tu- 
tukov, 1984). We show in Paper II that the time-scale 



t = 6.1 Gyr gives the correct mean color of normal galax- 
ies. For these galaxies, 36% of their total stellar mass 
formed below z = 1, 23% from z = 1 to 2, and 40% at 
z > 2. 

A feature dominant in the top panel is the strong spikes 
in the star formation history, indicating that the star for- 
mation takes place intermittently rather than continuously 
in normal galaxies. This feature is more manifest if we 
look at the star formation history of individual galaxies. 
Each galaxy experiences several conspicuous star forma- 
tion episodes, while they are quiescent for most of the time. 
This agrees with the suggestive evidence from the study of 
stellar ages in the Milky Way that it has expe rienced en- 



hanced episodes of star formatio n in the past (Majewsky. 



1993; Rocha-Pinto et al., 20001:). A study of star form- 



ing galax ies also supports the view t hat star formation is 
episodic ( Alonso-Hcrrcro ct al., 199(f ). 

The figure~[I ower panels) also shows that star formation 
ceased at high-redshift in small galaxies. Only 7% of stars 
form below z = 1 in galaxies with 2 x 10 8 < M s t e iiar < 
2 x 10 9 /i _1 M Q , and no stars form below z = 1 in galax- 
ies with Afstoiiar < 2 x 10 8 . This is consistent with the 
feature generic to the CDM structure formation, in which 
the formation of smaller structure ceases at higher redshift. 
Accretion of material stops, while most of residual gas is 
easily swept o ut by supcrnovae and photoion i zation from 
small systems (|Dckel fc Silk, 1986j; |Recs 1986|; |Efstathiou " 



1992|; IQuinn Katz, fc Efstathiou, 1996!; INavarro fc Stem 



mctz, 1997 ; Weinberg, Hcrnquist, & Katz, 1997: Gncdin 



2000b). Our prediction is that, in a CDM universe, smal 



galaxies found today are dominated by old stars; dwarf 
galaxies consist dominantly of stars of 10 Gyr old and do 
not show star formation activity in the recent epoch. It 
is known that most of dwarf galaxies in the Local Group 
are dominated by old ste llar population with occas i onal 
mixture of younger stars (Grebel, 1997 ; Mateo, 1998; [var] 
den Bcrgh, 1998: Gncdin, 2000a| ). A small admixture of 
young stars may not be excluded in the CDM model, since 



our calculation adopts a simplified treatment where only 
the evolution of the bulk of the material is followed, and it 
has no power to resolve minor components. The bulk star 
formation in small systems at late times, however, can- 
not be understood within the the scenario of hierarchical 
galaxy formation based on the CDM model. Major star 
formation activities in recent epochs, such as those in Leo 
I and Carina, do not fit the model 
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Fig. 2. — Star formation history of the three different samples of 
galaxies divided by their stellar mass at z = 0. Note the different 
scales on the ordinate. The dashed lines .are the least-square fits 
to the exponential function given in Table pj. Massive galaxies (top 
panel) continue to form stars until the present epoch, while less 
massive galaxies (bottom two panels) stop forming stars at higher 
redshift. 



5. MEAN AGE OF STARS AND GALAXIES 

The mean stellar age at each epoch in the history of 
the universe is also interesting. In Figure ||, we plot the 
median (solid square) and the mass-weighted mean age 
(open square) of all stars at each epoch. The median was 
calculated from the cumulative distribution of stellar mass 
at each epoch. The solid bars indicate the quartiles, which 
almost coincide with 1-er ranges shown with dotted bars. 
The average age of the extant stars decreases only slowly 
with increasing redshift. The mean stellar age today is 
4 — 5 Gyr and the age becomes smaller than 1 Gyr only 
at z = 3. The age 4 — 5 Gyr is slightly younger than the 
observed median stellar ag e of a 6 Gyr in the Milky Way 
( |Rocha-Pinto et al., 2000bj ). We show in Paper II that the 



distribution of B — V colors of simulated galaxies agrees 
with that of the local galaxy sample. The median age of 
1 Gyr at z — 3 means that the majority of galaxies at 
this redshift are dominated, in addition to star forming 
activities, by Balmer absorption features of A stars, which 

3 We note in passing that, in the warm dark matter scenario, low 
mass galaxies form primarily by "top-down" fragmentation rather 
than "bottom-up" accumulation 
in the standard CDM picture ([Bode et al., 2000 



arc found later than 
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would be tested with near infrared spectroscopy for Lyman 
break galaxies. 




Redshift 

Fig. 3. — Median (solid squares) and the mass-weighted mean 
stellar age (open squares) of all the stars at each epoch. The solid 
bars show the quartiles, and the dotted bars show The median 
and the quartiles were obtained from the cumulative distribution of 
the stellar mass as a function of the age at each epoch. The mean 
stellar age decreases slowly with increasing redshift, and exceeds 
1 Gyr at 2 rj 3. 

Figure f| gives the mean formation redshift of galaxies 
versus their stellar mass at z = 0. The mean formation 
redshift was calculated by taking the mass- weighted mean 
of the formation epoch (in time unit) of constituent stellar 
particles, and then converting it to redshift. Each dot rep- 
resents individual galaxy. The solid square symbols are the 
median, and the solid bars are the quartiles in each mass 
bin. This figure confirms the result we saw earlier that 
less massive galaxies found at z — are older than mas- 
sive ones. Small systems formed at high and intermediate 
redshift merge into more massive systems at lower red- 
shift. The mean formation redshift mostly decreases with 
increasing galaxy stellar mass, but it turns over slightly at 
the very massive end: this is due to the merging of oldest 
dwarfs into very massive galaxies. 

6. THE METALLICITY DISTRIBUTION AND THE 
AGE-METTALICITY RELATION 

The metallicity of stars as a function of their formation 
time (going from right to left) is plotted in Figure || for the 
same subsamples as in Figure ^|. Each point in the figure 
corresponds to one stellar particle in the simulation. The 
contours of an equally spaced logarithmic scale in number 
density is used where the density of data points is large. 
The plot shows that the stars which formed at low-redshift 
[z <^ 1) do not show unique metallicity but distribute from 
0.1-Z© to I.OZq. The width of the distribution slowly in- 
creases with redshift, while the upper envelope stays at 
1.0 Zq up to z rs 1. The widening trend of the age disper- 
sion is conspicuous at redshift z > 2. Note that the widen- 
ing takes place dominantly towards lower metallicity, while 



the upper envelope slightly increases to super-solar values. 
At z ;> 3, the distribution ranges from 10~ 6 Zq to 3.0Z Q . 
It is interesting to observe that the median of metallicity 
changes very little at z <J 2. The mass- weighted mean 
metallicity varies also slowly as 0.55,0.45,0.61,0.24, and 
0.22Z© for z = 0, 1, 2, 3, and 5, respectively. The decrease 
of mean metallicity above z = 2 is clear, whereas we can 
not conclude the change of mean metallicity for z < 2 
because of small statistics. 

It is known that stars in the Milky Way do not obey 
a unique age-metallicity relation, but show a significantly 
wide distribution ranging from O.IZq to 3Z Q . Whether 
the dispersion increas es with redshift, howeve r, is still 
a matter of deba te. Edvardsson et al. (1993, see also 
McWilliam (1997)) argued that the dispersion increases 
as stellar age gets older, while the metal rich stars form 
an envelope that is about solar metallicity independent o f 
stellar age. On the other hand, Rocha-Pinto ct al. (2000a) 
argued that age and metallicity show a stronger correla- 
tion, earlier stars being metal poor. Our result agrees with 
Edvardsson's, but not with Rocha-Pinto's result. 

The result we discussed above me ans that early galaxies 



are not ne cessarily metal deficient. Kobulnicky and Zarit 



sky (1999) show that emission line galaxies at z — 0.1 — 0.4 



have metallicity only slightly lower th an solar, which is not 
uncommon in local galaxy samples. Pettini et al. (2000) 
show that Lyman break galaxies at z w 3 have m 0.25Z Q , 
only modestly sub-solar. The strength of metal lines in the 
highest redshift quasar ( Fan et al. , 2000 ) is also suggestive 
of normal metallicity at z = 5.8. These observations are 
all consistent with our simulation. 

It is also our prediction that the Milky Way contains 
highly metal poor stars. We expect that 1% of stars have 
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Fig. 4. — Mean formation time (converted to redshift) of galax- 
ies vs. stellar mass of galaxies at z = in the simulation. Mean 
formation time is calculated by taking the mass-weighted average 
of formation time of constituent particles of each galaxy. The solid 
square symbols show the median and the solid bars show the quar- 
tiles in each mass bin. Less massive galaxies formed earlier, and 
later merged into massive ones. 
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Fig. 5. — Metallicity distribution of stars for the same subsamples 
as in Figure pi Each point in the figure corresponds to one stellar 
particle in the simulation with different masses. The abscissa is the 
age of the universe going from right to left, and the contours are 
of the number density distribution of stellar particles on an equally 
spaced logarithmic scale. .-The mass-weighted mean metallicity at 
each redshift is given in § 0. Stars formed earlier (z > 4) have wide 
range of metallicity ranging from 10 -6 — I.OZq, while those formed 
at z < 1 have roughly 0.1 — I.OZq. 



<J 10~ 4 Zq. Those stars are necessarily old. Another inter- 
esting result is that most super-metal rich stars are from 
high redshift epoch. This would be consistent with the 
fact that most of super- metal rich stars arc in the bulge 
rather than in the disk ([Frogel fc Whitford, 19871; |Rich." 



19880 if the bulge contains more old stars than the disk 
as usually conceived^. At the same time, our calculation 
predicts that the mean metallicity of the old population is 
lower than that of the younger population, which is also 



consistent with the observation (Edvardsson ct al., 1 993 
McWilliam & Rich, 1994j; |Rocha-Pinto ct al., 2000a|J 



We may ascribe the change of scatter in metallicity from 
high to low redshift to the following reason. Meta llicity is 
a strong fu nction of local overdensity , as shown by Gncdin 
and Ccn and Ostrikcr (1999c). Both authors have 



shown that the metallicity distribution in high-density re- 
gions is narrower than that in low-density regions. As 
time progresses, the universe becomes more clumpy, and 
star formation takes place only in moderate to high over- 
density regions where gas is already polluted by metals, 
resulting in a narrower scatter of metallicity at late times 
as seen in the top panel of Figure ^. 

Let us now consider the G-dwarf problem; the fact that 
there are too few metal poor stars in the solar neighbor 



hood compared to th e closed box model (van den Bergh 
1962|n |Schmidt, 1963[ ). In our model, a large fraction of 
stars formed in an early epoch, and such a population 
contains a significant fraction of metal poor stars. A his- 

4 A contradictory view is presented by Rocha-Pinto ct al. (2000g] ), 
who claim that super-metal rich stars are young 



togram analysis of the metallicity distribution at z — 
shows that 50% of stars have Z < 0.3Z@ in our simula- 
tion. This is roughly the same as the prediction of the 
closed box model, whereas the observation shows the frac- 



tion of these stars to be only a fe w percent ( Pagel fc Patch- _ 
ctt, 1975| ; Sommcr-Larscn, 1991). The qualitative trend of 



our prediction docs not seem to be easily modified. We 
ascribe, however, the presence of large metal poor stars 
in our simulation to our inability of resolving the internal 
structure of galaxies (i.e. disk from bulge). If the disk is 
a later addition to the bulge as the observations suggest 
(Fukugita, Hogan, & Peebles, 1996), it would contain lit- 
tle metal poor stars, as is clear from Figure || above. The 
material accreted onto galaxies at lo w-redshift was already 
pollu ted by metals from bulge stars ( |Ostrikcr and Thuan. 



1975). A high fraction of metal poor stars could occur in 



the CDM model, when we do not distinguish the disk and 
the bulg e components. Th is view is consistent with the 
work by Kauffmann (1996), where she provides a natural 
explanation to the G-dwarf problem, using a semianalytic 
model, by pre-enrichment by early outflows from dwarf 
progenitors. 

Finally in Figure |^, we show the mean metallicity of 
galaxies as a function of their stellar mass in the simula- 
tion. The mean metallicity of each galaxy is calculated 
by taking the mass-weighted average of constituent stel- 
lar particles. Solid square points are the median, and 
the solid bars are the quartiles in each mass bin. Mas- 
sive galaxies (M stc iiar > 10 10 /i _1 M©) at z = have 
mean metallicity of 0.1 — I.OZq, while less massive galax- 
ies have a wide range of mean metallicity ranging from 
10 -4 — 1.0-Z©. There are a few dwarf galaxies which have 
even lO -6 ^ outside of the figure. The star formation 
in less massive galaxies (identified at z = 0) took place 
with a wide range of metallicity at high-redshift as we saw 
in the bottom panel of Figure |5|, while larger galaxies at 
z = had star formation continuing up to the present 
epoch with normal metallicity of 0.01 — I.QZq. The two 
solid lines in the figure are obtained from a fit to the ob- 
servational d^a_of_sj3iral_jmd^^ given in 
Figure 4 of Kobulnicky and Zaritsky (1999). We converted 
the absolute £>-magnitudc to the stellar mass using the re- 
lation logMsteiiar = -0AM B + (2.75 ± 0.20), which the 
majority of galaxies in our simulation follow before dust 
extinction (see Paper II), and scaled the metallicity by 
\og(Z/Z e ) =12 + \og{0/H) - 8.89. A good agreement 
is seen between the simulated result and the observation 
in th e range of 10 7 < M stc iiar < lO 10 /?." 1 M Q . We remark 
that Kauffmann (1996] ) obtains a result similar to ours, 
which again suggests that the reason we are not seeing a 
similar cumulative distribution to her result and the so- 
lar neighborhood observation is simply because we cannot 
resolve the bulge and the disk into separate components. 

7. CONCLUSIONS 

We have studied predictions of the ACDM universe with 
regard to the star formation history and the stellar metal- 
licity distribution in galaxies. Our purpose has been to 
clarify what observations would test the validity of the 
model of galaxy formation based on ACDM. We have first 
shown that the global SFR averaged over all galaxies de- 
clines with a characteristic time-scale of r ss 6 Gyr. About 
a quarter of stars formed by z « 3.5 and another quarter 
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Fig. 6. — Average metallicity of galaxies vs. mean stellar mass 
at z = in the simulation. Each point represents one galaxy. The 
median in each mass bin is shown as a solid square, and the solid 
bars with it are the quartiles. The two solid lines in the figure are 
obtained from the best-fit to the obaecaationa] ^ata nf Qpiral anH j 



LII- 

regular galaxies given in Figure 4 of Kobulnicky and Zaritsky (1999|). 



See text for details. Massive galaxies (Ai stc n ar > 10 lu /i~ 1 Mq) have 
normal metallicity of 0.1 — I.OZq, while less massive galaxies take 
a wide range of values (10~ 4 — I.OZq). 



licity of dwarf galaxies today ranges widely (l(r 4 -lZ ), 
which reflects the large scatter of metallicity of stars 
formed at high-redshift. 

The G-dwarf problem in the solar neighborhood offers 
an interesting insight for formation of galactic structure. 
If we apply the prediction of the metallicity distribution 
to solar neighborhood, we encounter the G-dwarf problem: 
the prediction of the metallicity distribution does not differ 
much from that of the closed-box model. In our case, 
however, the G-dwarf problem would be solved if bulges 
are early entities and disks are later additions. 

We have briefly discussed observational tests for each 
item of the predictions. In many aspects the predictions 
are supported by observational evidence, but for some as- 
pects the predictions do not seem to agree with the current 
interpretation of the observations. The currently available 
observations, however, are sometimes controversial among 
authors or do not offer unique interpretations free from 
various working assumptions. We consider that decisive 
tests will be a task for the future. 

We thank Simon White for suggesting that we check the 
stability of the galaxies in our simulation. We arc grateful 
to the anonymous referee and Michael Strauss for useful 
comments on the draft. K.N. is supported in part by the 
Physics Department of Princeton University. M.F. is sup- 
ported in part by the Raymond and Beverly Sackler Fel- 
lowship in Princeton and Grant-in- Aid of the Ministry of 
Education of Japan. R.C. and J. P.O. are partly supported 
by grants AST 98-03137 and ASC 97-40300. 



by z « 2. This star formation history is consistent with 
the empirical Madau plot if modest dust obscuration is 
assumed. 

Our calculation indicates that star formation in L* 
galaxies continues intermittently to the present epoch as 
they accrete gas and merge with smaller systems with a de- 
cline rate corresponding to an e-folding time of t w 6 Gyr. 
Star formation in less massive galaxies ceases earlier with 
an e-folding time of r < 1 Gyr. In particular, dwarf galax- 
ies cease their star formation by z = 2. It is not excluded 
that minor star formation activity at low-redshift is missed 
in our calculation due to yet insufficient resolution, but it 
is unlikely that the global trend would be modified. We 
consider that this is a feature generic to the CDM structure 
formation scenario, and observational tests would provide 
a vital test. 

Stars formed in the CDM galaxy formation scenario do 
not follow a unique age-metallicity relation, but show a 
considerable spread of metallicity for a given age of stars. 
Metallicity of very young stars spreads from 0.1 — I.OZq. 
The spread gradually increases towards high redshift. It 
takes 0.01 - 1.0Z o at z = 2. At z = 3, the distribu- 
tion spreads over from 10~ 6 Zq to super-metal rich (3Zq). 
Young stars are necessarily metal rich, but old stars are 
not necessarily metal poor. The variation of the mean 
metallicity with cosmic time is only gradual. The average 
metallicity for galaxies at z « 1 differs little from that of 
local galaxies. 

On the other hand, mean metallicity varies as a function 
of the mass of galaxies. Less massive galaxies are metal 
poor. Another interesting prediction is that mean metal- 
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